Observation of quantum interference as a function of Berry's phase 
in a complex Hadamard optical network 
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Emerging models of quantum computation driven by multi-photon quantum interference, while 
not universal, may offer an exponential advantage over classical computers for certain problems. 
Implementing these circuits via geometric phase gates could mitigate requirements for error cor- 
rection to achieve fault tolerance while retaining their relative physical simplicity. We report an 
experiment in which a geometric phase is embedded in an optical network with no closed-loops, 
enabling quantum interference between two photons as a function of the phase. 



When a quantum mechanical system evolves under 
some Hamiltonian, the probability amplitudes associated 
with indistinguishable events can accumulate dynamical 
and geometric phases [l[ and interfere constructively or 
destructively. In Hong Ou Mandel interference Q , when 
two photons meet at the input ports of a beamsplitter, 
the event that both photons are transmitted is indis- 
tinguishable from the event that both photons are re- 
flected, but the associated probability amplitudes have 
opposite phases so interfere destructively: the probabil- 
ity to detect one photon at each output port is zero. This 
quintessentially quantum photonic interference generates 
the non-classical correlations in multi-photon quantum 
walks 0, 01 and the computational complexity of many- 
photon interference in large optical networks [5|-|7[ • These 
emerging models of quantum computation are unlikely 
to be universal, but may be exponentially more powerful 
than classical computers for certain problems. Crucially, 
since the basic models do not require initial entangle- 
ment, conditional gates, or feed-forward operations, large 
scale examples will be substantially less challenging to 
physically construct than a universal quantum computer. 
Achieving fault tolerance in these schemes without sacri- 
ficing their relative physical simplicity to unwieldy error 
correction is a key goal. 

Geometric phases and, more generally, non abelian 
holonomies have been proposed 0jS &- method to imple- 
ment fault-tolerant gates for universal quantum compu- 
tation 0-11 1, since they are robust against perturbations 
to which the important global geometric properties are 
invariant [l2l - [l7j . As described by Berry H, a phase is 
accrued by a vector in an instantaneous eigenstate of a 
Hamiltonian undergoing cyclic adiabatic evolution. An- 
ticipated in an earlier classical result [18[ and verified ex- 
perimentally the geometric phase has undergone 
important generalisations, including the non- abelian 22 1 



and non-adiabatic cases [23|, the non-cyclic [24l - l26j and 
non-unitary cases, and the case where the endpoints 
of the evolution are orthogonal psI - lBOj . The geometric 
phase has been observed at the single photon level and 
in the context of non- locality [31n33j , while a bi-photon 
wavepacket in a superposition of modes in a closed inter- 



ferometer exhibits the predicted increase in sensitivity 



to a geometric phase [3J-l36j that is observed for a dy- 
namical phase [37|. To date, however, all observations 
of optical geometric phases involve self-interference of a 
single photon or classical light interference. 

In light of the computational attributes of quantum 
interference between different photons and the desire 
to achieve fault tolerance in physically feasible com- 
putational models driven by this effect, demonstrating 
exquisite control over photonic quantum interference via 
an intrinsically robust geometric phase gate is a key step. 
Such an experimental connection in the context of these 
models is somewhat analogous to the implementation of a 
holonomic two-qubit gate in the established circuit model 
of universal quantum computation. Furthermore, to di- 
rectly observe the influence of the geometric phase on in- 
terference between photons, any measurement statistics 
should not be obfuscated by other phase dependent phe- 
nomena. In particular, single-photon interference, which 
has already been demonstrated to be predictably recep- 
tive to the geometric phase, should ideally be indepen- 
dent from the geometric phase in such an observation. 

Here, we establish an experimental functional relation- 
ship connecting a variable geometric phase (vGP) to si- 
nusoidal quantum interference between individual pho- 
tons of a pair. The vGP is imparted inside a four mode 
optical network that contains no closed loops, such that 
no single-photon interference can take place. Applied to 
only one photon of the pair in one of the modes, the vGP 
arises through a traversal of the polarisation-sphere com- 
prising a closed cycle and an open path; the end points 
of the total traversal are mutually orthogonal. The other 
three modes traverse lengths on polarisation-sphere equal 
to that of the vGP mode, but these include a partial, or 
total, path retracing such that a fixed GP, or no GP, is 
finally imparted. Polarisation vectors evolve via paral- 
lel transport, ensuring zero dynamical phase here, while 
the unknown dynamical phase contributions from small 
physical length mismatches are fixed. We observe high 
visibility quantum interference fringes, and find an ap- 
proximate flat line response for one-photon inputs, con- 
firming the absence of single-photon interference. 
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Complex Hadamard matrices [38| relate the computa- 
tional basis of a discrete Hilbert space to some mutually 
unbiased basis (soj . The four mode complex Hadamard 
unitary, H4, transforms quantum states according to 
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so that a single particle prepared in a well defined po- 
sition corresponding to an element of the computational 
basis, when acted upon by a device described by Hi, 
emerges with maximal uncertainty in its position. A 
large ensemble of similarly prepared one-particle input 
states will be found after the device with approximately 
one quarter of their total number at each of the four de- 
tectors (in the case of no losses). Modulation of the phase 
9 in the H4 device has no consequence for the maximal 
uncertainty in position of the single particle, and no con- 
sequence for the detection statistics of the ensemble. 

The optical network shown schematically in Fig. [TJa) 
consists of four one-half reflectivity beamsplitters, a swap 
of the two middle modes, a phase shift 9 on the lowest 
mode, and four detectors (-Dj); it is equivalent to H4 
with the labelling of input and output ports indicated. 
The network contains no closed-loop interferometers and 
a photon injected into any input port emerges in an equal 
superposition of the four output ports, regardless of the 
phase setting. Similarly, measurement in the computa- 
tional basis of one-photon ensembles cannot reveal any 
information about 9, and these statistics should ideally 
give a flat line response to modulation of this phase. 

The situation is dramatically different for two-photon 
input states. Simultaneous injection into the complex 
Hadamard network of a photon in mode |0) and a pho- 
ton in mode leads to a state that experiences pho- 
tonic quantum interference, producing correlations be- 
tween pairs of detectors as a function of 9. The condi- 
tional probabilities for coincidental detection of photons 
are summarised in Eq. ([1]). Given a detection at Di, the 
probability Pr(j | i) for a detection Dj is 

Pr(i ± 1 (mod 4) I i = {0,2}) = (1 ±cos6i)/2 

Pr(2 + 2|2 = {0,1}) = 0. (1) 

Implementing 9 geometrically creates a system in which 
quantum correlations from two-photon inputs are a func- 
tion of a geometric phase, whereas its effects are un- 
obscrvable in statistics from one-photon inputs, which 
would not be the case in an interferometer such as the 
Mach Zehnder. 

Devices implementing discrete instances of transforma- 
tions similar to have been constructed in bulk optics 
for a small number of dynamical phase settings (4ol | and 
in multimode waveguides (Zlj where the phase is per- 
manently fixed at a single value. Here, we encode in 
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FIG. 1: Complex Hadamard network and vGP. (a) Schematic 
of the optical network, with phase shift 6. (b) Experimental 
construction based on beam displacers, with waveplates im- 
plementing a geometric phase 9' . (c) Sphere showing stages 
of polarisation vector travel in the 9' section. 



position (rail) and polarisation, using a pair of paral- 
lelised Jamin-Lebedeff interferometers in a calcite beam 
displacer architecture, as shown in Fig.[ljb), to construct 
a device equivalent to H/i up to trivial phases on input 
and output ports, and relabelling (50| . Expanding to six 
modes through the 9' section, we use an arrangement of 
waveplates to implement the phase geometrically, allow- 
ing continuous transition between all phase values. The 
action of the beamsplitters in Fig. [IJa) is equivalently 
implemented on polarisation by half waveplates (A) with 
optic axes set to 22.5°; the swap is facilitated with beam 
displacers (B) that force horizontally polarised light to 
walk off at an angle, while vertically polarised light con- 
tinues undeviated; polarisation beamsplitters {S) convert 
polarisation to position for detection. 

The 9' section in Fig. [Ijb) comprises quarter wave- 
plates with optic axes fixed at 45° (Q), or at —45° (Q^), 
and half waveplates (P) free to rotate but with optic axes 
locked together at the same angle a. Overall, this imple- 
ments a polarisation flip on the top rail (t), the identity 
on the middle rail (m) |5l| , while the bottom rail (6) ex- 
periences a polarisation flip and a vGP as a function of a. 
Figure dKc) shows the polarisation-sphere [53- Photons 
in rails t and m, retrace their polarisation steps from the 
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FIG. 2: Photonic quantum interference and one-photon response to the change in experimental geometric phase 9' . (a) Raw 
data for all six coincidence outputs showing quantum interference as a function of AS', (b) The four AS' sensitive quantum 
signals for a 2-k range of the geometric phase, normalised, (c) Raw data for all eight output signals for both individual 
one-photon inputs. 



halfway point; while the total path length traversed on 
the polarisation-sphere for light in rail h is equal to that 
of rails t and m, the crucial difference in the actual route 
traversed leads to the experimental vGP, 9' {a). 

The polarisation vector in rail h, of light travelling the 
9' section of Fig. [ijb), makes a traversal of the sphere of 
Fig. [TJc) that includes a closed cycle and an open path. 
The full traversal is VRP(a)LHRP(a)LH, with P(a) in- 
dicating a point on the sphere determined by the variable 
angle of the half waveplates. The respective SU{2) uni- 
taries for the polarisation vectors of the three rails are 
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where \H) = |0) and \V) = |1) in each subspace. In this 
construction, a 1° change in the synchronised optic axes 
of both rotatable half waveplates results in a 4° shift in 
the experimental geometric phase: 4Aa = /S.9' . Light 
in rails t and h receive a polarisation flip during the full 
traversal so that the initial and final polarisation vectors 
of the respective rails are orthogonal. Interestingly, a 
geometric phase is accrued in this situation [28| and can 
be identified as the fixed i factor on these unitaries. 



To observe the correlations predicted in Eq. ([T]), pairs 
of photons were injected into the network as shown in 
Fig.HJb): a horizontally polarised photon into rail m, and 
a vertically polarised photon into rail 6, corresponding 
to computational states |0) and |1) respectively. The 
P waveplates were rotated almost two full revolutions, 
scanning a near 8 x 27r range for 9' . To confirm the 
insensitivity of one-photon statistics to /S.9' , one photon 
was injected into the network, for each of the |0) and |1) 
inputs, with another photon sent directly to a detector, 
as a herald. Photons were generated in a spontaneous 
parametric down conversion source [ssf . 

Experimental data shown in Fig. [2] strongly support 
predictions from Eq. ([T]) . Figure EJa) displays raw data 
for all six coincidence counts, with the expected four high 
visibility quantum interference fringes — the difference 
in amplitudes is due to different coupling efficiencies. In 
contrast, the two signals from detectors Dq & D2 on rail 
t (upper pane in Fig.[2{a)), and Di & D3 on rail m (lower 
pane in Fig.[2l^a)) show the predicted continuous and near 
total destructive quantum interference, with negligible 
response to A6''; the greatest phase response of these two 
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signals was Di k, D3 with a 0.9 ± 0.1% amplitude [54 1. 

Coincidence counts for the four A9' quantum response 
curves were re-taken over the 2tt range shown in Fig. ^h) 
and normalised for the measured coupling efficiencies. 
Following the notation from Eq. ([T]), the upper pane in 
Fig. shows P(l I 0) and P(l | 2) with a respective 
minimum and maximum at tt radians; the lower pane in 
Fig. [2lb) shows P(3 | 2) and P(3 | 0) with a respective 
minimum and maximum at tt radians. The average vis- 
ibility of these fringes is 94.0 ± 0.2% given from the fit 
that is shown as a solid black line. The small flat data 
in Fig. El^b) are respective accidentals signals, which are 
taken into account for the visibility calculation [55[ . The 
raw data in Figure EJc) confirm the insensitivity of one- 
photon statistics to A9'. The four panes on the left of 
Fig. ^c) show plots taken with a photon input in state 
0) with the top plot showing data for detector Dq, the 
next plot below showing data for Di and so on; similarly 
the plots on the right of Fig. ^c) are taken for input 
state In each case the average count is shown as the 
solid line and the phase insensitivity is quantified by the 
relative standard error (RSE) from this line. The average 
RSE from all eight plots is 3%; the plot with the greatest 
RSE at 6% is from detector D3 with input state 

The polarisation operations in the 9' section of 
Fig. [TJb) give zero dynamical phase contributions if the 
polarisation vectors are parallel transported at all times. 
The condition for parallel transport, for a certain state 
vector is {'ip{t)\tp{t)) = 0, which implies that, at 

an infinitesimal time step later, the sate \il}{t + St)) re- 
mains in phase with \'ijj(t)). This condition is met for in- 
dividual single photon states and for the composite two 
photon state superposed across all three rails of the 9' 
section in Fig. W[h). The lengths of rails t, m, and b are 
not matched on the wavelength-scale of the photons so 
that unknown dynamical phase differences occur between 
these rails, which are fixed due to the intrinsic stabil- 
ity of the experimental architecture. Therefore, the only 
variable phase is that on rail b and is a function of the 
common adjustable angle of the two P half waveplates. 

We have observed photonic quantum interference 
fringes that are a function of a variable geometric 
phase. This direct observation of the Berry phase in 
the quantum signal is possible because the optical com- 
plex Hadamard network contains no closed loops, so does 
not support single-photon interferometry. All active ele- 
ments of the state's Hilbert space make equal path length 
traversals on the polarisation sphere, but only those el- 
ements of the two-photon state that traverse via a par- 
ticular route, lead to a variable geometric phase. This 
route comprises a full cycle and an open path on the 
polarisation sphere, with mutually orthogonal start and 
end points. We have therefore simultaneously experimen- 
tally tested several important and distinct aspects of the 
geometric phase in conjunction with making the central 
quantum-interference observation of the Berry phase. 



The experimental circuit reported here is a small scale 
example of a model of quantum computation driven by 
photonic quantum interference 0, 01 , with a holonomic 
component. Generalisations of holonomies have previ- 
ously been widely investigated for robust gate operation 
in qubit based universal quantum computation, as dis- 
cussed, and theoretically examined in the specific case of 
photonic qubits Aside from imperfect unitary op- 

eration, the other major source of error in current linear 
optical experiments is typically photon loss, however, ev- 
idence suggests that loss in these models can be compen- 
sated through the injection of higher numbers of photons 
[4^. Furthermore, holonomies may reduce the exposure 
of photons to mode mismatch, leading to less demand on 
filtering, thus reducing loss. 

Scaling up examples of holonomic multi-photon- 
interference-drivcn computational models in waveguides 
(4^ where dynamical logic gates have been shown to work 
with high fidelity j45j is an appealing prospect. Any uni- 
tary transformation on modes can be implemented with 
a network of Mach Zehnder interferometers [i^ which 
have been realised in waveguides with variable thermo- 
optic 47 1 and electro-optic [i^ phase shifts, and inte- 
grated into a partially reconfigurable on-chip logic gate 
[49| . An interesting line of research is to consider the class 
of holonomic operations available given a large-scale fully 
reconfigurable optical unitary network, and the extent to 
which the global properties are invariant to imperfect 
splitting ratios in directional couplers and small random 
fluctuations from voltage controlled phase shifters, which 
act locally. 
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